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ABSTRACT 

The  gap-coupled  acoustoelectric  convolver  developed  at  Lincoln 
Laboratory  is  a surface-acoustic-wave  (SAW)  device  consisting 
of  a LiNbO-j  delay  line  and  a silicon  strip  supported  on  a series 
of  spacer  rails  (or  posts)  which  have  been  ion-beam  etched  into 
the  LiNbO-j  surface.  The  silicon/air-gap/LiNbCXj  structure  forms 
an  over-moded  acoustic  waveguide.  A theoretical  model  has  been 
developed  which  analyzes  the  perturbing  effect  of  rails  aad  pre- 
dicts the  mode  structure  and  beating  phenomena  between  modes. 
The  scattering  by  support  posts  is  also  analyzed. 
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MODAL  ANALYSIS  OF  SAW  CONVOLVER 

I.  INTRODUCTION 

1 2 

In  surface-acoustic-wave  convolvers,  ' a semiconductor  is  placed  in  close  proximity  to  the 

piezoelectric  medium  so  that  the  electric  field  of  the  SAW  can  interact  with  the  semiconductor 

to  produce  the  desired  nonlinear  effect.  Two  types  of  supporting  structure  have  been  tried  in- 
3 

eluding  the  random-post-support  and  the  finite-periodic  rail-support  structure.  The  post  sup- 

4 

port  is  unsatisfactory  because  of  the  excessive  loss  and  nonuniformity  found  among  the  devices 
tested.  The  rail  structure'’  is  a more  recent  development  which  seems  to  have  several  advan- 
tages over  the  random-post-support  structure.  Because  it  is  a simpler  and  more  deterministic 
structure,  the  effects  of  parameter  variations  are  more  controllable  and  suitable  for  analytical 
study. 

Indeed,  measurements  by  Reible*’  on  the  prototype  rail-support  device  show  some  very  in- 
teresting phenomena  such  as  the  laser-probe  scan  of  the  acoustic-field  profile  (Fig.  1)  and  the 
frequency  response  of  the  output  (Fig.  2).  A distinctive  feature  in  these  measurements  is  the 
beat  characteristic.  As  evident  from  Fig.  1,  the  acoustic  field  repeats  spatially  in  a periodic 
fashion.  We  will  refer  to  the  periodicity  as  the  beat  length.  Experimentally,  the  beat  length 
was  measured  to  be  approximately  137  Rayleigh  wavelengths.^  The  output  frequency  response 
of  Fig.  2 shows  a frequency  ripple  in  the  order  of  3 dB.  The  beat  is  the  result  of  interference 
between  the  acoustic  modes  of  the  device  which  are  analogous  to  the  modes  of  a periodic  loaded 
transmission  line.  The  discontinuity  presented  by  rail  structure  causes  the  incident  Rayleigh 
wave  to  be  scattered  periodically  in  the  transverse  direction.  This  scattering  process  resem- 
bles the  reflection  of  a wave  in  a transmission  line  by  discrete  elements  spaced  periodically  at 
fixed  intervals  along  the  line. 
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Fig.  2.  Output  frequency  response. 


Fig.  3.  Equivalent  circuits  in  a 
transmission  line:  (a)  two  port, 
(b)  series  element,  and  (c)  shunt 
element. 


In  adopting  the  periodic  transmission  line  model,  we  have  assumed  first  of  all  that  the  finite 
array  of  rails  can  be  approximated  by  an  infinite  periodic  structure.  This  seems  to  be  a good 
assumption  in  view  of  the  acoustic-field  profile  measurement  which  shows  no  changes  of  the  en- 
velope structure  connected  with  edge  effects  of  the  finite  transverse  dimension.  Secondly,  the 
finite  width  of  the  rails  is  ignored  by  modeling  them  as  discrete  elements  in  a transmission  line. 
This  approximation  is  justified  because  the  width  is  small  compared  with  the  transverse  Ray- 
leigh wavelength. 

In  Sec.  II,  we  review  the  scattering  matrix  of  two  ports  and  properties  of  periodic  struc- 
tures. Section  III  gives  a preview  of  the  kinds  of  modes  to  be  expected  in  the  periodic  rail 
structure  by  providing  a qualitative  description  of  the  guided  modes  and  their  dependence  on  the 
reflection  coefficient.  Then  in  Sec.  IV,  we  give  the  elastodynamics  of  the  rail  structure.  The 
analysis  of  this  section  enables  the  characterization  of  the  rail  structure  in  terms  of  its  scat- 
tering matrix  coefficients.  It  is  shown  that  the  derived  scattering  matrix  obeys  the  conserva- 
tion of  energy.  Using  the  model  developed  in  Sec.  IV,  we  next  study  in  Sec.  V the  reflection  co- 
efficients and  the  modes  in  the  periodic  rail  structure  as  a function  of  the  geometry  of  the  device. 
Finally,  in  Sec.  VI,  the  problem  of  Rayleigh  wave  scattering  by  the  post  structure  is  solved. 

II.  REVIEW  OF  PERIODIC  TRANSMISSION  LINE 

In  this  section,  we  will  briefly  review  the  scattering  matrix  description  of  two-port  junc- 
tions and  discuss  the  periodic  transmission-line  structure  made  up  of  a cascade  connection  of 
basic  unit  cells  consisting  of  a two  port  and  a section  of  transmission  line. 


A.  Scattering  Matrix 

In  general,  the  incident-  and  scattered-wave  amplitudes  of  a two-port  junction  shown  in 
7 

Fig.  3(a)  are  related  by 


(1) 


where  a^  and  a2  are  the  incident  waves,  b^  and  b 2 are  the  reflected  waves  of  ports  1 and  2,  and 
S>11'  S12,  S21>  and  S22  are  the  parameters  of  the  scattering  matrix  S.  For  our  problem,  we 
will  be  interested  only  in  a reciprocal  and  lossless  two  port  connected  to  transmission  lines  of 
equal  characteristic  impedance  at  both  terminals.  Reciprocity  and  symmetry  require  S^2  = S21 
and  = S22,  respectively;  losslessness  gives  the  conservation  of  energy  relationship 


lSll!2  + I Sl2l  2 = 1 • (Z) 

|S^|2  is  the  power  reflection  coefficient  and  |S^2|2  is  the  power  transmission  coefficient. 
Hence,  the  magnitude  of  S-2  is  uniquely  related  to  the  magnitude  of  S.^.  Furthermore,  it  can 

7 

be  shown  the  condition 


S11S12  + S12S22 


(3) 


always  holds  for  lossless,  reciprocal  two  ports. 


1 


Equation  (3)  relates  the  phase  of  S22  to  the  phase  of  S^. 


S 


11  = 


In  fact,  let 


then 


e2  = 91  + f T n,r  ‘ (4) 

To  summarize,  for  the  problem  we  are  investigating,  the  scattering  matrix  description  of 
a two-port  junction  consists  of  only  two  independent  parameters:  namely,  the  magnitude  and 
phase  of  the  reflection  coefficient  S^. 

The  basic  ideas  of  scattering  are  illustrated  by  considering  two  simple  examples.  In 
Fig.  3(b),  a series  reactance  (jX)  (normalized  to  the  line  impedance)  is  connected  in  a transmis- 
sion line.  The  incident  wave  is  labeled  a^  and  the  reflected  wave  as  bt  = S'  (assume  a2  = 0). 

We  can  think  of  the  transmitted  wave  (b2 ) as  made  up  of  two  parts:  a^  which  exists  in  the  ab- 
sence of  scattering  by  the  series  element  and  the  scattered  wave  which  arises  because  of  scat- 
tering. The  scattered  waves  (S  ) and  (S  ) can  be  thought  as  generated  by  an  equivalent  voltage 
source  (V)  induced  by  the  incident  wave.  To  see  how  the  equivalent  voltage  source  is  related 
to  the  scattered  wave,  we  write  the  following  equations: 

al  + bl  " b2  = <a!  - bi>  jX  (5a) 

or 

S'  - S+  = V = (&1  - b^  jX  (5b) 

and 

al  " bl  = b2  <M 


or 


-S'  = S+  . (6b) 

Equation  (5b)  indicates  that  the  voltage  source  is  the  difference  of  the  scattered  wave  on 
either  side  of  the  series  element.  Furthermore,  according  to  Eq.  (6b),  the  scattered  wave  on 
one  side  is  the  negative  of  the  other  side.  Combining  Eqs.(6a)  and  (5a)  results  in  the  familiar  ex- 
pression for  the  reflection  coefficient: 


s - bl  - ix 

S11  " ' T+lX 


(7) 


The  transmission  coefficient  is 


; - h - a + Si  - 2 

12  " a4  at  ' 2 + jX 


(8) 


The  1 in  the  expression  for  S^2  is  due  to  the  complete  transmission  of  the  incident  wave  a ^ in 
the  absence  of  scattering  by  the  series  element.  Obviously,  conservation  of  energy  [Eq.  (2))  is 
satisfied. 
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A similar  analysis  for  the  shunt  element  with  normalized  susceptance  jB  in  a transmission 
line  of  Fig.  3(c)  gives 


S + ST  = -I  (9) 

where  I is  the  equivalent  current  source.  Note  in  the  shunt  connection  the  current  source  is  the 
sum  of  the  scattered  wave  with 

S = S . (10) 

The  reflection  coefficient  is 

S - h - ~iB 

^11  " a4  " 2 + jB  ' (11) 

transmission  coefficient  is 


b2  S+ 
S,,  = -r-  = 1 + — = 


iZ~  al  ~ at  - 2 + jB  ' (12) 

and  conservation  of  energy  is  satisfied.  We  will  encounter  an  analogous  problem  in  the  scatter- 
ing of  a Rayleigh  wave  by  the  rail  structure. 


Fig.  4.  Periodic  transmission  line. 


B.  Periodic  Structure 


Having  reviewed  the  scattering  matrix  representation  of  two  ports,  we  next  take  up  trans- 
mission line  periodically  loaded  with  such  two  ports  as  shown  in  Fig.  4. 

The  wave  amplitudes  (C's)  are  related  by 


n m 


where  p is  the  propagation  constant  and  the  A matrix  is  defined  by  the  equation  which  relates 
the  wave  amplitudes  (C's): 

r +T  + i 

C A,,  A,,  C\. 

n li  12  n+1 

= (14) 

C"  A,.  A,,  C" 

nj  L 21  22J  n+lj 

and  A is  related  to  the  scattering  matrix  § by 


A - — 
A11  ' S„ 


A21  " A12 


S2  -S2 
b12  bll 


According  ic  Floquet's  theorem,  the  wave  amplitude  at  the  (n  + 1)  ST  terminal  can  differ  from 
the  wave  amplitude  at  n**1  terminal  by  at  most  a phase  delay.  Thus,  we  assume 

«„■«  ■ ^ 

where  ud  is  the  phase  delay.  Substitute  Eq.  (16)  into  (13)  gives  an  eigenvalue  equation  for  k. 

cos  *d  = j~-  [e37Jd  + (S22  - S24)  e"^d]  . (17 

12 

The  characteristic  modes  in  a particular  unit  cell  of  the  periodic  structure  is  given  by 


, , . if  ^n*  / 1 ^nd  S12  j«J\  ^nXl  f 

5n  x>  = 2 le  + (s^  e " S^e  ) 6 ] f°r 


°<x<| 


1 -^n(x'd)  . (J_  ^nd  fl2  j..d\  ^n(x-d) 

= 2 l ls22  e e ) e 


j/3n(x-d)i 


for  < x < d . (18) 

Solution  of  Eq.  (17)  for  a typical  periodic  structure  gives  the  familiar  dispersion  diagram 
as  shown  in  Fig.  5.  It  shows  "bands"  of  alternating  even  and  odd  symmetry  as  exhibited  by  the 
mode  pattern  at  the  Brillouin  zone  center  (k  = 0)  and  edge  ( k = v/d).  These  two  points  are  of 
particular  interest  to  us  because  we  accept  as  modes  only  those  solutions  which  have  vanishing 
group  velocity.  With  a uniform  excitation  of  the  periodic  rail  structure,  we  expect  from  sym- 
metry that  the  acoustic  wave  is  the  superposition  of  standing  waves  in  the  transverse  direction. 
Furthermore,  because  of  the  assumed  uniform  transducer  input,  only  symmetric  modes  such 
as  modes  1 and  5 could  be  excited.  We  will  refer  to  these  modes  as  the  fundamental  mode  and 
the  higher-order  mode.  These  two  modes  are  most  strongly  excited  and  the  beating  between 
them  gives  rise  to  the  observed  spatial  and  frequency  ripples. 


© 

® 

(2)  © I EVEN 


Fig.  5.  Dispersion  diagram  and  modes 
of  periodic  structure. 


TRANSVERSE  DIRECTION 

III.  QUALITATIVE  DESCRIPTION  OF  THE  MODES 

Before  going  into  the  detailed  modeling  of  the  rail  structure,  we  now  discuss  qualitatively 
the  behavior  of  the  fundamental  and  higher-order  modes  as  a function  of  reflection  coefficient. 

The  fundamental  and  higher-order  mode  patterns  are  shown  in  Figs.  6(a-d)  as  the  magnitude 
and  phase  of  the  reflection  coefficient  are  varied.  The  magnitude  of  S^  is  reduced  by  50  per- 
cent in  Fig.  6(b)  as  compared  with  that  of  Fig.  6(a).  As  expected,  weaker  reflection  enables  the 
wave  function  <p  (x)  to  spread  out  from  the  space  between  the  reflecting  elements.  The  same  ef- 
fect can  be  achieved  by  retarding  the  phase  of  the  reflection  coefficient  as  evident  in  Fig.  6(c) 
in  which  phase  is  decreased  by  40  percent  from  that  of  Fig.  6(a).  A dramatic  spreading  effect 
results  as  shown  in  Fig.  6(d)  if  both  magnitude  and  phase  are  decreased. 

The  spreading  effect  makes  the  fundamental  mode  more  uniform  and  therefore  matches 
better  the  uniform  input  transducer  profile.  Since  the  higher-order  mode  is  orthogonal  to  the 
fundamental,  it  becomes  less  excited.  Therefore,  we  can  infer  immediately  a very  important 
fact:  in  a rail  structure  which  has  a reflection  coefficient  with  small  magnitude  and  retardation 
in  phase  at  the  mode  angles  corresponding  to  the  fundamental  and  higher-order  mode,  the  spa- 
tial and  frequency  ripple  size  would  be  small. 

Over  the  range  of  50-percent  variation  in  magnitude  and  40-percent  variation  in  phase,  the 
mode  patterns  basically  retain  their  shapes  (i.e.,  nearly  one  half-cycle  of  sinusoid  for  the  fun- 
damental and  three  half-cycles  for  the  higher-order  mode).  This  indicates  that  the  mode 
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pattern  is  not  strongly  dependent  on  the  reflection  coefficient.  Its  shape  is  inherently  a charac- 
teristic of  the  periodicity. 

For  the  range  of  variation  in  reflection  coefficient  studied,  the  beat  length  varies  from  80 
to  154  Apj  assuming  a periodicity  of  9 Xjj.  As  the  modes  spread  out  more  uniformly,  the  trans- 
verse wave  vector  becomes  smaller,  forcing  the  longitudinal  wave  vectors  of  the  modes  to  ap- 
proach more  closely  with  each  other,  therefore,  the  beat  length  tends  to  increase. 

Now  suppose  the  spacing  is  reduced,  but  assume  the  reflection  coefficients  remain  unchanged 
for  the  fundamental  and  higher-order  mode.  Then  the  transverse  wave  vectors  increase  lin- 
early, resulting  in  shorter  longitudinal  wave  vectors;  but  the  length  difference  is  greater,  hence, 
a smaller  beat  length.  In  reality,  reducing  the  spacing  means  increasing  the  mode  angle  and, 
thus,  the  reflection  coefficients  do  not  remain  constant  but  in  fact  would  decrease  causing  the 
modes  to  spread  more  uniformly.  It  is  not  clear,  without  a more  detailed  model  of  the  ridge 
structure,  which  is  the  dominant  effect  on  the  beat  length.  As  shown  in  Sec.  V,  the  dependence 
of  reflection  coefficient  on  mode  angle  is  not  strong  enough  to  overcome  the  dominant  longitu- 
dinal wave  vector  length-difference  effect.  The  net  result  of  decreasing  the  spacing  is  to  de- 
crease the  beat  length.  At  the  same  time,  since  the  modes  are  more  uniformly  spread  out,  the 
ripple  size  would  decrease. 

IV.  SCATTERING  MATRIX  REPRESENTATION  OF  THE  RAIL  STRUCTURE 

With  the  framework  of  periodically  loaded  transmission  set,  we  are  now  ready  to  study  in 
more  detail  the  physics  of  the  rail  structure.  In  particular,  we  are  interested  in  characterizing 
the  scattering  of  a Rayleigh  wave  by  me  rail  structure  (Fig.  7)  in  terms  of  its  scattering  matrix. 
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The  basic  idea  in  solving  the  scattering  problem  is  to  find  an  equivalent  source  of  scatter- 
ing, in  our  case  a stress  tensor,  in  place  of  the  line  discontinuity.  The  stress  tensor  must  sat- 
isfy boundary  conditions  specified  by  the  dynamics  of  the  rail  in  the  presence  of  the  incident 
Rayleigh  wave.  For  example,  suppose  the  rail  is  pinned  so  that  no  motion  is  allowed.  This 
can  be  realized  by  assuming  an  infinitely  massive  rail  or  an  infinitely  stiff  superstrate.  Then, 
to  satisfy  the  pinned  boundary  condition,  there  must  be  an  equivalent  stress  at  the  position  of 
the  rail  causing  a displacement  which  exactly  cancels  the  displacement  of  the  incident  wave  and 


l 


of  the  scattered  wave  as  sensed  by  the  source.  The  effect  of  the  scattered  Rayleigh  field  must 
be  included  in  the  figuring  of  equivalent  sources  because  the  scattered  field  reacts  back  on  the 
source  which  produced  it  in  the  first  place.  This  feedback  phenomenon  is  encountered  in  other 
scattering  problems  as  well  and  is  very  important  if  conservation  of  energy  is  to  be  satisfied. 
Notice  that  we  have  allowed  only  the  surface  wave  in  this  feedback  process.  The  bulk  wave  is 
ignored  because  it  is  not  found  experimentally  to  be  significant/1 

To  find  the  source  which  satisfies  the  near-field  boundary  condition,  we  must  know  for  a 
given  point  stress  what  is  the  displacement  in  the  neighborhood  of  the  source.  In  other  words, 
we  like  to  evaluate  the  dynamic  Green's  function  near  the  origin  where  the  source  is  located. 

As  it  turned  out,  the  dynamic  Green's  function  within  a distance  from  the  origin,  which  is  small 
compared  with  the  wavelength,  is  identical  to  the  static  Green's  function.  An  analogous  situa- 
tion occurs  in  electromagnetic  problems.  The  near  field  of  an  oscillating  dipole  is  equal  to  the 
quasi-static  dipole  field.  We  will  justify  this  equivalence  in  more  detail  later  when  we  study 
the  Fourier  integral  expression  of  the  Green's  function.  The  asymptotic  far-field  evaluation  of 
the  Fourier  integral  at  the  Rayleigh  pole  gives  the  scattered  surface  wave. 
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Fig.  8.  Rail-support  structure, 
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A.  Green’s  Function  of  a Point  Source 


Suppose  stress  with  Fourier  transform  cr^fk)  is  applied  at  the  origin  of  the  x-y  plane  as 
shown  in  Fig.  8;  we  like  to  find  the  scattered  surface  wave  on  an  isotropic  medium  with  Lame 
constants  \,p.  Assume  the  potentials  (p  and  ty,  respectively,  for  the  compressional  and  shear 
waves  of  the  form: 


where 


<p  = C e 
0 = jS  e 


;j(kxx+kyy)  „pz 


-i(k  x4k  v) 

J x yJ  ’ -sz 

J pt 


X 

, 2 


(19) 

(20) 


(21) 

(22) 
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(23) 


k2  = k 2 + k 2 

x y 


ki  = u)(r^2ii)1/2 


(24) 


and 


k2  * 

where  p is  the  mass  density. 

£ 

Suppose  a surface  stress  exists,  Fourier  decomposition  gives 


(25) 


C(k)  = 


(k  + s^)  ^zz(k)  + j2ks  [ CJxz(k)  cose  + cTyz(k)  sine] 
flD(k) 


(26) 


S(k) 


2kP<?zz(k)  + j(k2  + s2)  [<rxz(k)  cose  + a (k)  sine] 

per!] 


(27) 


where 


D(k)  = (k2  + s2)2  — 4k2ps  (28) 

and  e is  the  angle  the  line  of  observation  (y-axis)  makes  with  the  x-axis.  The  angular  depen- 
dence comes  from  the  projection  of  the  shear  stress  components  onto  the  line  of  observation. 
The  reason  is  that  only  the  shear  stress  pointing  along  the  direction  of  observation  couples  with 
the  Rayleigh  wave  in  that  direction.  Taking  superposition  of  all  spatial  harmonics,  we  have: 


<p(x,  y,  z)  = i dk  \ dk  C(k)  e j(kxx+kyy)  e~Pz 

- oo  — oo  y 

r°°  C°°  -j(k  x+k  v)  _ 

ip(x,y, z)  = j dkxJ  dky  iS(k)  e r e s 


(29) 

(30) 


For  a point  stress  with  the  total  applied  force  components  F , F , and  F 

x y £. 


zz=-Fz6(x)  6(y) 

(31) 

xz=-Fx6(x,6(y) 

(32) 

yZ  =-Fy  6(x)6(y) 

(33) 

the  G reen1  s function  is 


<?G(x,  y,  z)  = — !-*•  \ dk  1 
1 2n)  x J- 


[Zr)^  d _oo  *■  j.oc  y 

x.-i'VV’.-P* 


« - (k2  + s2)  F,  - j2ks(F  cose  + F sin©) 

dk  - i 


pD(k) 


itG(x,y,  z)  = — j-y  1 dk  l 

(Ztt)  ^ -oo  J -a 


dk 


— 2kpF  — j(k2  4 s2)  (F  cose  4 F sine) 


(2? r)1"  ^-oo  x ^-oo  y 

-j(kXX4kyy)  -S 


■pTW 


x e 


(34) 


(35) 
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Now,  using  Eqs.  (34)  and  (35),  we  justify  our  earlier  claim  that  the  dynamic  Green's  function 

near  the  source  appears  quasi-static.  Within  distance  of  the  order  of  a wavelength,  integration 

in  k-space  does  not  show  any  special  dependence  until  and  ky  are  much  greater  than  2tt/X ^ or 

larger  than  either  k.  or  k,.  In  this  range  the  integrand  in  Eqs.  (34)  and  (35)  is  independent  of  fre- 
1 q ^ G 

quency.  Therefore,  <p  and  ip  are  identical  with  the  static  Green's  function  in  the  neighborhood 
of  the  source. 

B.  Static  Green's  Function  of  a Point  Source 

The  static  Green's  function  for  the  displacement  due  to  a point  stress  is  given  by  Landau 

9 

and  Lifshitz. 


Ux  = TTE  7 [ --  r-*-  Fz  + 2(1  - Fx  + 7T  (xFx  + yFy>] 

Uy  = 1171  7 | — ~r'?)y  Fz  + 2(1  - Fy  + ^ (xFx  + ^1 

| [2(1  - a)  Fz  + (1  - 2a)  I (xFx  + yFy)| 


u = 1 

z 2nE 


(36) 

(37) 

(38) 


where  E is  the  Young's  modulus,  a is  the  Poisson  ratio,  and  r is  the  distance  from  the  source. 
Equations  (36)  to  (38)  pertain  to  a point  stress  and  therefore  are  not  directly  applicable  to  our 
rail  problem  which  has  a line  geometry.  However,  as  we  will  see  later,  these  equations  are 
directly  applicable  to  the  post  problem.  In  the  next  section  we  will  derive  the  counterparts  of 
Eqs.  (36)  to  (38)  for  the  rail  problem. 

C.  Green's  Function  of  a Modulated  Line  Source 

Figure  9(a)  illustrates  the  geometry  of  scattering  for  the  rail  structure.  The  rail  is  as- 
sumed to  lie  along  the  y-axis  with  the  Rayleigh  wave  incident  at  an  angle  0 with  respect  to  the 
x-axis  and  the  reflected  wave  going  off  also  at  angle  0. 

We  will  assume  an  equivalent  modulated  line  source  along  the  y-axis.  The  modulation 

-j'V  -3kvy 

e J comes  from  the  fact  that  the  incident  Rayleigh  wave  has  e J dependence.  Analogous 
to  Eqs.  (34)  and  (3  5),  we  take  a Fourier  superposition  of  space  harmonics  now  only  in  the  kx 
dimension. 


G,  . 1 r~  „ — (k2  + s2)  F — j2ks(Fx  cos0  + F sin©)  -j(k  xtk^l 

V (x.y.s)  = 2i  j_oo  dkx  jn5?Ej e e 

,y,z)  = ^J  c 

*->  — oo 


( (x 


-2kpFz  — j(k  + s ) (F  cos0+  F sin0)  -j(k  x+k  v) 

^d?E5  1 e e 


(39) 


(40) 


where  the  modulated  line  source  is 


yz 


-ikvy 

e y 6(x) 

(41) 

-jkyy 

e 6 (x ) 

(42) 

e y 6(x) 

(43) 

13 


The  integration  in  Eqs.  (39)  and  (40)  can  be  done  easily,  by  picking  up  the  Rayleigh  pole  residue 
at  kx  = J kR  - ky  where  kR  = 2tt/xr  is  the  propagation  constant  of  the  Rayleigh  wave;  the  inte- 
grations are  done  in  Appendix  A.  The  results  are  repeated  here  for  convenience. 


G,  v -2kRs(Fy  sine  ±Fx  cose)  +j(2k^-k22)Fz  -jkRr  _pz 
v {r’z>  = jTD^^Tcose e e 


, G,  -J(2kR  -k2)  (F  sine*  Fx  cose) -(2kRP)  Fz  -jkRr_gz 


(r,  z)  = 


pD'(kR)  cos e 


The  Green's  functions  for  displacement  Ur  and  Uz  of  the  Rayleigh  wave  on  the  surface  (z  = 0) 


are  then 


tt  _ dtp  _ dip 

r “ ?r  3z 


= ^<kR)  cose  [kR(2kR  - k2  ~ 2Ps)  Fz  + 3k2  s(Fy  sin0  * Fx  COS0)I 


U = ^ +11 
uz  3z  + 3r 


^(kR)  cose  fjk22PFz  - kR(2kR  ~ k22  - 2PS)  (Fy  sin0  * Fx  cos0)l 


To  simplify  notation,  Eqs.  (46)  and  (47)  are  written  in  tensor  form: 


U.  = E..F.  e 

i 13  3 


2 

±jk2  s cose 
E11  = " pD'(kR) 


2 

jk2  s sine 


-12  - pD'(kR) 


kR(2kR~k22-2Ps) 


pD'(kR) 


2 

± jk2  s sine 


-21  - pD'(kR) 


, ..  2 . 2„ 

+]k2  s sin  e 


-22  - pD' (kR)  cose 


E„  = -E 


23  ' 32 


2 2 

^2  sin0 


nD'(k^)  cos 0 


l 


I 

j 

j 


' 


*kR(2k^  - k22  - 2ps) 
pD'(kR) 


"■33  " pD'(kR)  cos© 


D.  Stiffness  Matrix 

The  near  field  of  a modulated  line  source  approaches  that  of  the  quasi-static  Green's  func- 
tion. We  assume  the  force  is  uniformly  distributed  across  the  width  of  the  rail  (a  good  assump- 

-ikvy 

tion  if  the  width  « A R)  and  sinusoidally  modulated  with  a phase  factor  e J along  the  length  of 
the  rail.  The  displacements  at  the  center  x = 0 are  found  by  convolutions  of  Eqs.  (36)  to  (38). 
This  is  done  in  Appendix  B.  The  result  is  given  in  the  form 

Wj  l57) 

where  T^  is  a diagonal  tensor  which  measures  the  stiffness  of  the  medium. 

E.  Reaction  of  Scattered  Wave  on  the  Source 

As  in  any  acoustic  scattering  process,  the  scattered  wave  tends  to  react  back  on  the  source 
of  scattering.  The  interaction  of  the  scattered  wave  with  the  source  was  illustrated  in  Sec.  II 
using  the  examples  of  series  and  shunt  elements  in  a transmission  line.  Now  we  want  to  see 
explicitly  how  this  interaction  occurs  in  the  case  of  scattering  of  a Rayleigh  wave  by  the  rail 
structure. 

A great  deal  can  be  deduced  from  the  symmetry  of  the  problem.  For  example,  Fx  excites 

only  a surface  wave  of  odd  symmetry  with  respect  to  the  y-axis,  while  Fz  and  F excite  surface 

waves  of  even  symmetry.  Therefore,  Fx  does  not  "feel"  the  scattered  wave  generated  by  Fz  or 

F and  vice  versa.  Stated  in  another  way,  F does  not  couple  with  F or  F ; however,  F and 
y x z y z 

Fy  do  couple  with  each  other.  This  observation  suggests  the  scattering  of  a Rayleigh  wave  can 
be  viewed  as  being  made  up  of  an  antisymmetric  scattering  in  the  x-direction  and  a symmetric 
scattering  in  the  y-  and  z -direction.  Indeed,  as  we  shall  see  later,  this  turns  out  to  be  the  case. 

Let  denote  the  average  of  the  scattered  wave  which  reacts  back  on  the  source.  Then 

we  can  relate  <(U^>  to  F^  by  the  following  matrix  equation. 
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0 

0 

0 

V22 

V23 

0 

V32 

V33 

VU=Ell(+)  • 


V22 = E22 


V23  = V32  = E23 


V33 = E33  * 
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The  "averaging"  is  taken  so  that  Fx  "feels"  the  difference  of  the  scattered  wave  (analogous  to 
Eq.  (5b)  for  scattering  by  a series  element  in  a transmission  line),  while  F and  F "feel"  the 

y z 

sum  of  the  scattered  wave  [analogous  to  Eq.  (9)  for  scattering  by  a shunt  element  in  a transmis- 
sion line). 


F.  Dynamic  Equation  of  the  Rail  Structure 

The  boundary  condition  satisfied  by  the  sources  F^  obey  the  dynamic  equations: 


U: 


.K  + T.PfP  + <U.>  = T.Pf.’ 
i D 1 N r i]  3 


U^U 


(59) 


and 


-(F.L  + F,P)  = — 2ahp  u)2U. 
J J S J 


(60) 


Referring  to  Fig.  7,  T,.  and  T..  are  the  stiffness  tensors  of  the  upper  and  lower  substrate 
l3  L3  Tj  t 

with  applied  force  per  unit  length  F.  and  F.  , respectively,  2a  is  the  width,  h is  the  height. 


R - 


3 


p is  the  mass  density  of  the  spacer,  U.  is  the  incident  Rayleigh  displacement,  and  U . is  the 

“ Ip  o J 

net  displacement.  The  displacements  13  ^ at  the  surface  z = 0 are  given  by 


K - 


-3(kvx+k„y) 


2kr 


x '~y  cose 


(61) 


R k2  -3<k  x+kyy) 

Uy--iA+2i^e  ^ sine 


(62) 


R 2kR  ~k2  ~ 2sp  -3<kxx+kyy) 


Uz=A+  2s 


(63) 


Substitution  of  Eqs.  (60)  and  (58)  into  Eq.  (59)  yields 


<Tij  + T£  - 2ahps*2TG<TF  + Vlj)  + V^J  F,  = -(^  - 2ahps<u2TJ)  uf 


(64) 


where  6y  is  the  Kronecker  delta.  Let 


Gh  = + - 2a6p  2T,P (T,P  + V,  J + V„ 


'I  13 


13 


il ' lj 


(65) 


and 


R..  = -(a..  — 2ahp  o;2T.P) 
13  13  s 13 


(66) 


Equation  (64)  becomes 


^4 

O 

0 

0 ‘ 

"f 

X 

’Ru 

UR1 

X 

0 

G22 

G23 

Fy 

= 

R22 

0 

G32 

G33. 

.FZ. 

_R33 

(67) 


The  form  of  the  matrix  equation  (67)  confirms  the  fact  that  antisymmetric  scattering  in  the 
x-direction  is  independent  of  the  symmetric  scattering  in  the  z-  or  y-directions. 

G.  Scattering  Matrix  of  the  Rail  Structure 

To  find  the  scattering  matrix  of  the  rail  structure,  we  consider  an  incident  Rayleigh  wave  R 

D 

with  displacements  U.  launched  at  angle  0 with  respect  to  the  x-axis  against  the  rail  [see 
3 ST 

Fig.  9(a)].  The  reflected  wave  S and  transmitted  wave  T have  displacements  IE  and  tk  , re- 
spectively. Because  the  displacements  of  the  Rayleigh  wave  are  related  to  each  by  definite  ra- 
tios [see  Eqs.  (61)  to  (63)],  we  can  define  a set  of  unique  scattering  matrix  coefficients  from  the 
ratios  of  the  z-displacements  of  S and  T waves  to  the  z-displacement  of  the  R wave,  i.e.. 


U«,  U and  have  one-to-one  correspondence  with  voltages  a^,  b1>  and  b2>  respectively, 
of  the  transmission  line  model  shown  in  Fig.  3(a).  The  excitation  of  Fig.  9(a)  can  be  decomposed 
into  its  antisymmetric  part  [Fig.  9(b)]  and  symmetric  part  [Fig.  9(c)].  Indeed,  by  adding  the  dis- 
placement components  of  the  antisymmetric  excitation  to  that  of  the  symmetric  excitation,  we 
get  exactly  the  displacements  of  the  original  excitation  represented  in  Fig.  9(a). 

In  the  antisymmetric  excitation,  all  the  incident  Rayleigh  displacement  components  cancel 
except  the  x-component  which  is  equal  to  the  x-displacement  component  of  the  R wave.  Solving 
for  Fx  in  Eq.  (67),  substituting  into  Eq.  (48),  we  find  the  antisymmetric  scattered  field  normal- 
ized to  the  incident  field  to  be 

z 


1 4 jxd 


k.fs  (1  — 2ahp  (o^Tj^ ) cose 

X = pD'(k  ) L ,,  , _U  ,_L  77  2,„  U , 

R Tll<1  +Tll/Tll  ~2ahPs“  T11) 
Therefore,  we  find  for  the  antisymmetric  excitation: 


i 


The  — 1/2  in  Eq.  (72)  is  due  to  the  total  transmission  of  the  incident  wave  — R/2  from  the  bottom 
side  in  the  absence  of  scattering;  similar  reason  explains  the  l/2  in  Eq.  (73). 

For  symmetric  excitation,  the  x-component  of  the  incident  Rayleigh  wave  displacement  can- 
cels leaving  the  y-component  and  z-component  which  are  equal  to  the  y-  and  z-displacement 
components  of  the  R wave.  The  symmetric  scattered  field  is  found  by  first  solving  for  and 
Fz  in  Eq.  (67)  and  then  substituting  into  Eq.  (48),  we  find 

E32Fy  + E33Fz 


_ 


2 + jB- 


pD'(k_)  cose 

r\ 


2 

+ s sin  0 


Pip)- 


(1  +T2U2/T2L2  -2ahps<.2T3L3) 


(1  + T33/T33  -2ahpsa)2T3E  ) 


(1  +T2L2/T2L2  -2ah^sw2T2U2>  +T3VT3L3-2ahPs“2T2U2) 


Thus,  for  the  symmetric  excitation: 


TJ  W / o 

rs  = Hl.  , 1 + Ds  = 1 g.-3B_s 

U R 2 1 2 + jBs 


U ' ' s 

s uz  1 , „s  1 2 - iBs 

T = 5-  a 7 + D = ^ - — -- 

UcR  2 7 2 + jBs 


Finally,  the  scattering  matrix  coefficients  of  the  original  excitation  [Fig.  9(a)]  are  found  by 
superposition: 


Sll=  uR 
z 


u(sa)+  u(sS) 

Z Z 


and 


u (Ta>  + u (tS) 

z z 


= . (79) 

1 + jXa  2 + jBs 

H.  Conservation  of  Energy 

An  important  test  of  the  validity  of  and  given  by  Eqs.  (78)  and  (79)  is  to  show  that  the 
conservation  of  energy  is  satisfied. 


I jX*  _ ]BS  I2 
li+jXa  2+jBs| 


jXa  jBs  I2 

1 + jXa  2 + jBS> 


= i Q.E.D.  (80) 

This  condition  becomes  obvious  by  noting  T's  and  r's  of  Eqs.  (72),  (73),  (76),  and  (77)  have  mag- 
nitude of  l/2  implying  conservation  of  energy  holds  for  the  antisymmetric  part  and  symmetric 

V.  RESULTS  ON  THE  RAIL  STRUCTURE 

In  this  section,  we  present  calculations  of  the  reflection  coefficient  and  modes  of  the  peri- 
odic rail  structure  using  the  model  described  in  the  previous  section. 

A.  Reflection  Coefficient 

Figure  10  plots  the  fraction  of  reflected  power,  ( ^ | 2,  of  the  ridge  vs  the  angle  of  inci- 
dence 9.  As  expected,  | ^ | 2 increases  with  angle  of  incidence  and  approaches  unity  at  grazing 
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Fig.  10.  Reflection  of  Rayleigh  wave  by  the  rail  structure. 
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leu lated  from  the  assumed  Lame  constants, 


incidence  (e  = tt/2).  The  reason  is  because  the  ridge  appears  to  be  "softer"  for  a Rayleigh  wave 
incident  at  a small  angle  as  evident  by  looking  at  the  stiffness  matrix  t-  of  Eqs.  (B-10)  to  (B-12). 

The  mass  of  the  spacer  makes  little  difference  in  the  reflection  coefficient  as  shown  by  the 
curve  labeled  (2)  corresponding  to  the  case  of  LiNbO^,  Si,  and  zero  mass.  On  the  other  hand, 
the  scattering  of  the  Rayleigh  wave  is  very  sensitive  to  the  stiffness  of  the  superstrate.  For 
example,  by  letting  the  superstrate  become  infinitely  stiff  (curve  1 ) the  reflection  increases. 

On  the  other  hand,  removal  of  the  superstrate  (curve  5)  leads  to  a dramatic  decrease  in  the  re- 
flectivity. The  effect  of  the  superstrate  is  to  deaden  the  surface  vibration.  This  effect  is  en- 
hanced by  a wider  ridge  as  evident  in  curves  2,  3,  and  4 where  the  width  is  changed  from 
6 to  4 pm  and  then  to  2 pm.  Although  not  shown  in  Fig.  10,  the  phase  of  the  reflection  coef- 
ficient generally  decreases  with  a decrease  in  the  magnitude  except  for  the  case  of  no 
superstrate. 

B.  Modes  of  the  Periodic  Rail  Structure 

Using  the  scattering  matrix  coefficients  of  the  ridge,  we  proceed  to  find  the  modes  of  the 
periodic  rail  structure  using  Eq.  (17).  The  results  are  summarized  in  Table  I.  We  have  stud- 
ied the  behavior  of  the  fundamental  mode  and  the  higher-order  mode  by  changing  the  width  and 
the  spacing  of  the  rail.  The  resulting  mode  profiles  are  given  in  Figs.  11  to  14.  Beat  length 
and  frequency  ripple  are  also  calculated. 

The  mode  expansion  of  the  uniform  input  is  given  by 


0O 


2j  aiN.<p.(x)  = 1 
i=l 

(81) 

where  <p.(x)  is  the  mode  profile  given  by  Eq.  (18)  with  unity  peak  amplitude, 
tion  constant  such  that 

fU  is  the  normaliza- 

4 

(‘d  2 2 

l N i>“(x)  dx  = 1 

Jo 

(82) 

Ni=  \wt  sin/Ed)]-1/2 

(83) 

where  p.d  is  the  normalized  transverse  wave  number,  and 
pd  2N.  fS  d 

ai  = J N^fx)  dx  = sin  — . (84) 


The  modes  of  a 6.5-pm-wide  periodic  rail  structure  with  100-pm  spacing  are  shown  in 
Fig.  11.  Superposition  of  the  fundamental  and  higher-order  mode  according  to  Eq.  (81)  gives 
the  field  profiles  of  Fig.  12  depending  on  the  relative  phase  of  the  modes,  given  by  the  longitu- 


we  can  get  either  the  sum  or  the  difference  of  the  two  modes  resulting  in  a beat  pattern.  Spa- 
tially, at  a fixed  frequency,  the  superposed  field  profile  repeats  within  a distance  L^,  the  beat 
length,  such  that 

- 1 5)  = 2tr  . <86) 

The  beat  pattern  also  occurs  at  the  output  transducer  as  the  frequency  is  swept.  The  frequency 
interval  Af  is  related  to  by 


Af  _ ^b  (87) 

f " L 

where  L is  the  distance  between  the  input  and  output  transducer.  The  depth  of  ripple  as  de- 
tected in  the  received  power  as  a function  of  frequency  is 

, 2 2 

/al  ~ aS 
M = 20  log  I — 2 2 

\a!  + a5 

Table  I shows  that  for  the  6.5-pm  device,  the  beat  length  is  86  \R  as  compared  with  137  XR 
measured  experimentally6  (Fig.  1).  The  ripple  frequency  depth  is  -1.9  dB,  somewhat  less  than 
the  -3  dB  observed  in  the  output  frequency  response  (Fig.  2).  The  discrepancy  is  probably 
caused  by  groove  scattering  not  included  in  our  present  model. 

Figure  13  shows  the  modes  in  the  rail  structure  with  narrower  width  but  with  the  same 
spacing.  Because  the  reflectivity  is  reduced  and  phase  retarded,  the  modes  become  more  spread 
out  resulting  in  longer  beat  length  and  smaller  ripple  (see  Table  I)  as  explained  in  Sec.  III.  To 
take  advantage  of  the  mode-spreading  effect  with  weaker  reflectivity,  the  width  is  made  narrower 
and  the  spacing  is  reduced.  The  latter  allows  the  mode  angles  to  be  smaller  (i.e.,  closer  to 
normal  incidence),  therefore  achieving  a weaker  reflection  coefficient.  The  spacing  is  decreased 
for  another  important  reason.  By  increasing  the  number  of  supports,  the  load  supported  by 
each  rail  is  decreased;  as  a result,  the  deformation  is  kept  small,  minimizing  groove  scatter 
ing.  The  combined  effect  is  seen  by  the  mode  profiles  shown  in  Fig.  14  for  a structure  with 
2 -jim  width  and  30-^m  spacing.  It  is  clear  that  the  modes  have  indeed  spread  out  significantly 
as  compared  with  those  of  Fig.  11.  Consequently,  as  shown  in  Table  I,  the  depth  of  ripple  has 
reduced  to  -0.73  dB.  However,  Lb  has  shortened  to  11  XR. 

A word  of  caution,  by  continuing  the  reduction  of  the  rail  width,  eventually  groove  scatter- 
ing caused  by  larger  deformation  near  the  rail  becomes  important.  The  effect  of  groove  scat- 
tering is  certainly  to  increase  the  magnitude  of  reflection  coefficient.  However,  it  is  not  clear 
without  further  study  which  way  the  phase  will  change.  If  the  phase  is  retarded,  then  the  modes 
tend  to  be  spread  out  offsetting  the  effect  of  the  increased  reflectivity.  In  addition,  bulk  scat- 
tering may  now  play  an  important  role  as  the  groove  is  steepened  and  more  rails  are  in  place 
as  a result  of  reducing  the  spacing.  A definitive  judgment  cannot  be  made  until  we  study  more 
closely  the  coherent  bulk  scattering  by  an  array  of  rails. 

VI.  POST  SCATTERING 

In  this  section,  we  take  up  the  post  scattering  problem.  It  is  solved  in  a similar  fashion 
as  the  rail  scattering  problem.  We  first  of  all  evaluate  asymptotically  the  far  field  from  the 
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integral  expression  of  the  Green's  function  of  a point  source  derived  in  Sec.  IV.  Then,  the  stiff- 
ness matrix  and  the  dynamic  equation  are  derived.  Finally,  we  present  the  radiation  patterns 
of  post  scattering. 


J 


A.  Green's  Function 

Integral  expressions  of  the  Green's  function  are  given  in  Eqs.  (34)  and  (35).  In  Appendix  C, 
the  scattered  Rayleigh  wave  is  evaluated.  The  results  are 


<pG(r,  z)  = 


— 2kDs(F  sinG  + F cosO) 
R .y  x 

pD'(kR)/k,. 


+ j(2k2  -k22) 


x e 


-jkRr 


e’Pz 
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W 
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The  displacements  of  the  Rayleigh  wave  at  z = 0 are  then 

3^  _ 3^ 
r 3r  3z 
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(91) 


x [Fy  sine  + Fx  cose]  exp  [- j(kRr  - -| )] 


(92) 


B.  Stiffness  Matrix 

Assuming  that  the  forces  are  uniformly  distributed  over  a disk  of  radius  (b),  the  displace- 
ments at  the  center  are  found  by  convolution.  The  result  is  given  in  terms  of  the  stiffness 
matrix  T.j  (see  Appendix  D)  relating  displacements  to  the  forces  in  Eq.  (57). 

C.  Radiation  Patterns  of  the  Post  Structure 

The  dynamic  equation  for  the  post  structure  is  identical  in  form  to  the  rail  problem 
[Eq.  (59)]  except  that  we  ignore  the  reaction  of  the  scattered  wave  back  on  the  source.  The  dif- 
ficulty in  treating  the  reaction  of  the  scattered  wave  back  on  the  source  is  caused  by  the  fact 
that  a detailed  energy  balance  as  we  have  in  the  rail  problem  is  not  as  easily  carried  out.  First 
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of  all,  the  incident  wave  is  a uniform  plane  wave  of  infinite  extent,  while  the  wave  scattered  by 
the  point  source  is  more  or  less  radial.  Such  scattering  geometry  presents  difficulty  in  the 
normalization  of  energy  which  is  not  the  case  in  the  rail  scattering  problem.  In  addition,  the 
effect  of  the  scattered  wave  on  the  point  source  is  not  as  easily  treated  analytically  as  for  a 
line  source  because  of  the  nature  of  the  point  singularity.  ' Therefore,  it  must  be  kept  in  mind 
that  our  result  is  an  approximation  which  provides  only  a qualitative  picture  of  the  scattering 
process. 

Assume  a uniform  Rayleigh  wave  incident  on  the  post  in  the  x-direct ion.  The  induced  forces 
which  satisfy  the  dynamic  equation  are  given  by 


F = 
x 


Fz- 


-7Tb2(l  -p^T^) 
-7Tb2(l-pghu>2T3U3) 
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T - (1  + T3L3/T  - - Psh^T3u3) 


2rr.  U , 
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F = 0 
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(93) 


(94) 


(9  5) 


Substitution  into  Eqs.  (91)  and  (92)  gives  the  scattered  Rayleigh  wave. 

Note  the  scattered  field  decreases  as  1 /\fr  from  the  source  expected  for  two-dimensional 
scattering.  Furthermore,  the  radiation  pattern  of  the  vertical  source  component  Fz  is  isotropic 
while  it  has  a cosine  dependence  for  the  horizontal  source  component  F^.  The  displacement  ra- 
diation patterns  are  shown  separately  in  Fig.  15  for  a post  with  a silicon  superstrate.  It  is  inter- 
esting to  note  that  the  scattering  in  the  forward  direction  is  stronger  than  that  in  the  backward 
direction  because  of  the  constructive  interference  of  the  scattered  wave  components.  Further- 
more, we  observe  that  the  vertically  scattered  component  is  larger  than  the  horizontal  compo- 
nent. This  is  because,  at  the  surface,  the  Rayleigh  wave  displacement  in  LiNbOj  is  1.5  times 
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Fig.  15.  Displacements  of  scattered 
SAW  by  the  post  structure. 
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SUPERSTRATE 


Fig.  16.  Radiation  pattern  of  scattered 
SAW  by  the  post  structure. 


0 
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larger  in  the  vertical  direction  than  in  the  horizontal  direction,  inducing  a larger  F than  F . 

The  difference  is  further  enhanced  in  the  excitation  of  the  scattered  Rayleigh  wave. 

The  power  radiation  patterns  for  a post,  with  and  without  superstrate,  are  shown  in  Fig.  16. 
We  note  that  the  scattering  is  dramatically  increased  by  the  pinning  effect  of  the  superstrates, 
a phenomenon  also  observed  in  the  rail  scattering. 

VII.  SUMMARY 

We  have  presented  a model  for  the  scattering  of  the  Rayleigh  wave  by  the  rail-  and  post- 
support structure.  Using  the  scattering  matrix  coefficients  predicted  by  our  model,  we  have 
analyzed  the  modes  in  the  SAW  convolver  with  a periodic  rail  structure.  The  theoretical  results 
obtained  are  in  reasonable  agreement  with  the  experimental  measurements.0 

Our  study  provides  insight  into  the  optimal  design  of  the  rail  support  in  order  to  minimize 
the  moding  problem.  We  recommend  that  the  rails  be  made  thinner  and  placed  closer  together 
so  as  to  decrease  the  ripple  size  in  the  output  frequency  response. 

Finally,  we  recommend  further  study  into  the  scattering  of  the  Rayleigh  wave  by  rail-  and 
post-support  structure  by  taking  the  finite  dimension  of  the  support  structure  into  account.  This 
would  involve  a more  detailed  study  of  the  groove  scattering  as  well  as  the  bulk  wave  scattering 
problem. 

The  groove  scattering  may  very  well  account  for  the  discrepancy  between  the  theoretically 
predicted  frequency  ripple  using  our  present  model  and  the  experimental  measurement.^  Fur- 
thermore, in  order  to  optimize  the  rail -structure  design  and  to  get  a handle  on  the  pressure 
dependence  of  the  mode  pattern,  the  groove  scattering  problem  must  be  studied  carefully. 

Our  study  has  not  included  loss  due  to  the  bulk  scattering  which  is  expected  to  be  significant 
when  the  deformation  caused  by  the  rail  loading  becomes  larger.  Furthermore,  mr!e  conver- 
sion as  a result  of  scattering  by  nonuniform  loading  along  the  rail  and  scattering  by  boundary 
roughness  can  also  contribute  to  loss.  We  believe  our  model  is  still  incomplete  without  proper 
account  for  these  loss  mechanisms  which  are  important  in  understanding  the  limitations  on  the 
convolver. 
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APPENDIX  A 

GREEN'S  FUNCTION  OF  A MODULATED  LINE  SOURCE 
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where 


D(k)  = (k2  + s2)2  — 4k2ps 


(A-3) 


We  want  to  evaluate  integrals  in  Eqs.(B-i)  and  (B-2) 
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by  contour  integration  in  the  kx  complex  plane.  The  only 


2 2 

pole  singularity  occurs  at  kv  = ± J kR  - kv  . Letting  kR  be 


‘x  ^ ' R y ' 
slightly  imaginary  displaces  the  poles  off  the  real  axis  as 

shown. 

For  x > 0,  we  choose  to  close  the  contour  in  the  lower 
half -plane.  To  find  the  residue,  we  need  the  derivative  of 
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D(k)  evaluated  at  kx  = JkR  — k . This  is  done  by  expand- 


ing D(k)  in  the  neighborhood  of  kx: 
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where  / k 2 + k 2 + kR  and  R(k)  are  functions  with  weak  dependence  on  k as  k — /kR  - k 2 ; so 
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we  can  replace  the  two  functions  by  2kR  and  R(kR),  respectively,  continue. 
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Finally,  for  x > 0,  let  r = J x2  + y2  to  be  the  distance  measured  along  the  y-axis  (see  Fig.  8). 
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For  x < 0,  we  close  the  contour  in  the  upper  half -plane,  the  result  is  a change  of  sign  of 
the  Px  term‘  A similar  contour  integration  is  carried  out  for  i /iG.  In  summary, 
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where  the  plus  sign  is  used  for  x > 0 and  the  minus  sign  for  x < 0. 


APPENDIX  B 

STIFFNESS  TENSOR  FOR  THE  RAIL  STRUCTURE 


Convolutions  are  of  the  form: 


u.(x,  y)  = j dy*  j dx'  ^ u.9(x  - x'  , y - yr 


where  XJ/f  are  Green’s  functions  for  a point  source  given  by  Eqs.  (36)  to  (38)  and  2a  is  the  width 
of  the  rail.  We  are  only  interested  in  the  displacement  at  the  center,  i.e.,  x = 0.  Integrations 
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of  terms  with  U/T  odd  in  x'  vanish.  For  the  remaining  terms  we  need  the  following  integrals: 


r"  i -3kvy' 

■•'Lp  » dy' 

h-r 


I = r e‘3V'  dyI 

J -00  P 

I = r k-zpl  e"]  V dy, 

J -00  n 


Let  p = y'  - y 


P = vx'  + (y  - y' ) 


‘.~Vf.T= 


J-°°  / ,2  2 

Vx11  + p 

"■'kyy  ■»  r"  coskyK 

Jo  / ,2  2 

vx'  + p 


j-yr- 

e J dp 


dp  . 


Let  X = -4 
x' 


-jk__y  ,•>“  cos  (k  | 
It  = e r 2 ^ y 


0 VI  +r 


= e ^ 2KQ(k^|x’ |)  (Kq  is  the  modified  Bessel  function) 


For  kyX*  « 1, 


/ k |x’|  v -jky 
Ij  a - 2 I In  -J-j — 4 {I  e ^ 


[Eq.  (B-2)] 


where  f = 0.57721. 
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Let  x1  — 0,  we  note  I3  is  well  behaved  while  all  other  integrals  become  singular,  there- 
fore I3  can  be  ignored.  Next,  the  x' -integrations  are  done: 
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Now,  Eq.  (B-l)  can  be  put  in  the  form 
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where  the  stiffness  matrix  Ty  is  diagonal  with 
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APPENDIX  C 

GREEN'S  FUNCTION  OF  A POINT  SOURCE 
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Evaluating  the  kx  integration  as  done  in  Appendix  A gives 
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Steepest  decent  path  (S.D.)  is 
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